arXiv:hep-ph/9808378v5 17 Dec 1998 


Bremsstrahlung of a Quark 
Propagating through a Nucleus 


Boris Z. Kopeliovich^’^, Andreas Schafer^ 
and Alexander V. Tarasov^’^ 

^ Max-Planck Institut fiir Kernphysik, Postfach 103980, 69029 Heidelberg 
^ Joint Institute for Nuclear Research, Dubna, 141980 Moscow Region 
^Institut fiir Theoretische Physik der Universitat, Philosophenweg 19, 69120 Heidelberg 
'^Institut fiir Theoretische Physik, Universitat Regensburg, 93040 Regensburg 


Abstract 

The density of gluons produced in the central rapidity region of a heavy ion 
collision is poorly known. We investigate the influence of the effects of quantum 
coherence on the transverse momentum distribution of photons and gluons radiated 
by a quark propagating through nuclear matter. We describe the case that the 
radiation time substantially exceeds the nuclear radius (the relevant case for RHIC 
and LHC energies), which is different from what is known as Landau-Pomeranchuk- 
Migdal effect corresponding to an infinite medium. We find suppression of the 
radiation spectrum at small transverse photon/gluon momentum kx, but enhance¬ 
ment for kx > 1 GeV. Any nuclear effects vanish for kx > 10 GeV. Our results allow 
also to calculate the kx dependent nuclear effects in prompt photon, light and heavy 
(Drell-Yan) dilepton and hadron production. 

PACS numbers: 12.38.Bx, 12.38.Aw, 24.85.+p, 25.75.-q 
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1 Introduction 


One of the major theoretical problems in relativistic heavy ion physics is the reliable calcn- 
lation of gluon bremsstrahlung in the central rapidity region. It is one of the determining 
factors for the general dynamics of heavy-ion collisions, the approach to thermodynamic 
equilibrium and the possible formation of a quark-gluon plasma-like state. This problem 
has been approached by a variety of ways. We do not want to discuss the relative draw¬ 
backs and merits of the various approaches here and we will only cite those, which are 
directly related to ours. 

In this paper we consider bremsstrahlung of photons and gluons resulting from the 
interaction of a projectile quark with a nucleus for the case that the radiation time is 
much longer than the time needed to cross the nucleus. This radiation or formation time 
was introduced in [Q and can be presented as, 

cosh y 2a; 

fr = _ - ^ _ 

kx kx ’ 

where y, uj and kx are the rapidity, energy and the transverse momentum of the radi¬ 
ated quantum in the nuclear rest frame. Eq. m assumes that the radiated energy is 
relatively small, i.e. uj Eq. It is easy to interpret the formation time (|l|) as lifetime 
of a photon(gluon)-quark fluctuation p[ or as the time needed to distinguish a radiated 
quantum from the static field of the quark |^. 

The total time for bremsstrahlung is proportional to the initial energy and can there¬ 
fore substantially exceed the time of interaction with the target 0. Radiation continues 
even after the quark leaves the target. This part of radiation does not resolve multiple 
scattering processes. Important is only the total momentum transfer. This illuminating 
manifestation of coherence is along these lines that the well known Landau-Pomeranchuk- 
Migdal effect (LPM) for long formation times can be treated. Note that LPM effect corre¬ 
sponds to the opposite energy limit, when the radiation time is much shorter that the time 
of propagation through the medium. It was first suggested by Landau and Pomeranchuk 
m and investigated by Migdal ^ and has attracted much attention during recent years 
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i,B0, ii- This regime applies only for the problem of energy loss in a medium, which 
is not the problem we discuss here. Our treatment should apply to the real situation 
in heavy-ion collisions at high energies. The relationships between the cited papers are 
complex. In a recent publication Baier et ah 0 have shown that their diagrammatic 
approach is in fact equivalent to that of Zakharov [^]. The latter is, however, physically 
far more intuitive and therefore lends itself more easily to a generalization to the case that 


the nuclei are not inhnitely extended. In another recent paper Kovchegov and Mueller [11 


have undertaken the hrst attempt to calculate in-medium modihcation of the transverse 
momentum distribution of gluon radiation. This paper has also elucidated the relation 
between the approaches of [0 and 0. In the approach of [^ based on the use of the 
light-cone gauge the hnal state interactions summed up in [Q (in the covariant gauge) are 
effectively included in the light-cone wave function. These observations suggest that all 
three different approaches might be equivalent when followed carefully enough. 

The main goal of this paper is to study the dependence of the effects of coherence 
on the transverse momentum of the radiated photon or gluon. We use the light-cone 
approach for radiation hrst suggested in and developed in [^, As it is based on 
an explicit treatment of the transverse coordinates it is easily adapted to our purpose. In 
addition it seems to be by far the most direct and elegant approach. We described this 
approach in Section 2 for both photon and gluon bremsstrahlung. We establish a relation 
between the strength of the coherence effects and the transverse size of the Fock state 
containing the radiated quantum. 

The second main result of our paper is the extension of the light-cone approach to 
calculations for differential cross sections as functions of the transverse photon/gluon 
momentum kx- This is presented in section 3. As one might have expected, nuclear 
shadowing, i.e. suppression of radiation, is most pronounced at small kr- An unexpected 
result is antishadowing, i.e. enhancement of radiation for k^ > 1 GeV, which, however, 
vanishes for still larger 

The results and practical implications for the Drell-Yan process, prompt photon pro- 
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duction and hadropro duct ion are discussed in the last section. 


2 Integrated radiation spectra 


We start with electromagnetic radiation. We cover both, virtual photon radiation (dilep¬ 
ton production) and real photon radiation (so called prompt photons). 

The total radiation cross section for (virtual) photons, as calculated from the diagrams 
shown in Fig. |l|, has the following factorized form in impact parameter representation 


a 




Figure 1: Feynman graphs for bremsstrahlung. 


(see also 


0 ), 


\^^q{a,rT)f aq^iarr) . ( 2 ) 

Here T.yg(a, Ft) is the wave function of the y — q fluctuation of the projectile quark which 
depends on a, the relative fraction of the quark momentum carried by the photon, and 
tt, the transverse separation between 7 and q (T is not normalized). <Jqq{p) is the total 
interaction cross section for a qq pair with transverse separation p and a nucleon. crqq{p) 
depends also parametrically on the total collision energy squared s, a dependence we do 
not write out explicitly (see, however, section 4). This becomes only important when hts 
to actual data are performed. Eq.(2) contains a remarkable observation which is crucial 
for this whole approach [j^: although we regard only a single projectile quark, the elastic 
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amplitude of which is divergent, the radiation cross section is equal to the total cross 
section of a qq pair, which is hnite. 

This can be interpreted as follows. One should discriminate between the total interac¬ 
tion cross section and the freeing (radiation) cross section of a fluctuation. The projectile 
quark is represented in the light-cone approach as a sum of different Fock components. 
If each of them interacts with the target with the same amplitude the coherence between 
the components is not disturbed, i.e. no bremsstrahlung is generated. Therefore, the 
production amplitude of a new state (a new combination of the Fock components) is pro¬ 
portional to the difference between the elastic amplitudes of different fluctuations. Thus 
the universal divergent part of the elastic amplitudes cancels and the radiation amplitude 
is hnite. 

It is also easy to understand why the qq separation in (Q) is arx- As is pointed 
out above one should take the difference between the amplitudes for a quark-photon 
huctuation and a single quark. The impact parameters of these quarks are different. 
Indeed, the impact parameter of the projectile quark serves as the center of gravity for 
the 7 —g huctuation in the transverse plane. The distance to the quark in the quark-gluon 
Fock-state is then ar^p and that to the photon is (1 — a)r'p. 

The wave function of the 'j*q huctuation in (^) for transversely and longitudinally 
polarized photons reads (compare with [ 0 ), 

T^4(rT, a) = oTi Ko{erT) (3) 

2 71 

Here Xij cire the spinors of the initial and hnal quarks. Kq{x) is the modihed Bessel 


function. 

The operators 0^’^ have the form. 



0^ = i mqO^ e* ■ {n X a) + a e* ■ {a x V) — i{2 — a) e* ■ W , 

( 4 ) 


0^ = 2m^* (1 — a) , 

( 5 ) 

where 




-l- (1 — a)m^. . 

(6) 
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e is the polarization vector of the photon, n is a nnit vector along the projectile momentnm, 
and V acts on ty- For radiation of prompt photons m^* = 0. 

Eq. (H) can be nsed for nnclear targets as well. We consider hereafter formation times 
given by the energy denominator, 


tf — 


2 Eg a{l — a) 

+ rrig 


> R/ 


(7) 


which snbstantially exceed the nnclear radius. In this limit the transverse 'y* — q separation 
in the fluctuation is ’’frozen”, i.e. does not change during propagation through the nucleus. 
The recipe for the extension of Eq. (|^) to a nuclear target is quite simple [^, . One 

should just replace a^JarT) by a^Jarx), 


da^{q jq) 
(i(lna) 


= 2 y cPb J (Rtt \'^^q{a,fT)f |l — exp 


--aqg{arT)T{b) 


( 8 ) 


where 

/ OO 

dzpA{b,z) . (9) 

-OO 

Here pA{b,z) is the nuclear density which depends on the impact parameter b and the 
longitudinal coordinate One can eikonalize Eq. because a fluctuation with a ’’frozen” 
transverse size is an eigenstate of interaction [|^]. 

Eq. shows how the interference effects work versus kx- At small the exponent 
o'qq{<yrT)T{b)/2 -C 1 since aqq{ar'r) is small. Therefore, one can expand the exponential 
and the cross section turns out to be proportional to A. This is the Bethe-Heitler limit for 
bremsstrahlung. In the opposite limit aqq{arT)T{b)/2 3> 1 one can neglect the exponential 
for b < Ra and the cross section is proportional to . This is the limit of full 
coherence when the whole row of nucleons with the same impact parameter acts like a 
single nucleon. As the gluon transverse momentum is related to the inverse of r^, one 
could expect that the limit of maximal coherence is reached for small k^i and the Bethe- 
Heitler limit for large The situation is, however, more complicated as discussed in the 
next section. 


Gluon radiation is described by the diagrams [jT^ shown in Fig. The radiation cross 
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Figure 2: Feynman graphs for gluon bremsstrahlung of an interacting 
quark. 


section for a nucleon target and the nuclear effects [|I3 look similar to those of Eqs 


da^{q gq) 
d(\na) 


= 2 y d% J d'^VT |'hgg(a, — exp 


-7:(^gqqi.F,r2)T{b) 


( 10 ) 


where 'Fgg{a,rT) is the wave function of a quark-gluon fluctuation which has the same 
form as in Eq.(^), but with the replacements 7 * ^ g, aem 4as/3 and ^ mg. We 
keep the gluon mass nonzero in order to simulate the possible effects of conhnement on 
gluon bremsstrahlung. agqq is the interaction cross section of a colorless gqq system with 
a nucleon [p!7[| , 


(J. 


gqq 


9 r 'I 1 

(Fl, F2) = -|cXg-q(ri) + aqg{r2)\—- aqq{n - F2) , 


( 11 ) 


where Fi and F 2 are the transverse separations gluon - quark and gluon ~ antiquark 
respectively. In the case of gluon radiation, i.e. Eq. (0). Fi = Ft and F 2 = (1 — a)rT. 


Although Eq. (W) looks simple, it includes the effects of quark and gluon rescattering 
in the nucleus to all orders. 


3 The transverse momentum distribution 
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3.1 Electromagnetic radiation 


The transverse momentum distribution of photon bremsstrahlung in quark-nucleon inter¬ 
actions integrated over the hnal quark transverse momentum reads (see Appendix A), 

(Pa^{q gq) 1 


(i(lna;) (Pkr 
where 




exp[!A^T(f'i - 'f' 2 )l n) fi) tr,(ri, fq, a) , 


CT 


( 12 ) 

( 13 ) 


f(fi, f2, a) = - <Tqq[a{ fi - rq) 

By integrating over kr one obviously recovers Eq. (||), since (j^(r, r, a) = aqq{ar). 

For a -C 1 one can use the dipole approximation for the cross section, i.e. one can set 
<Jqq{p) = Cp^. Moreover, this approximation works also rather well at larger interquark 
separations, even for hadronic sizes ii. For the latter the cross section is proportional 
to the mean radius squared. Therefore, we use the dipole approximation for all cases 
considered. Then (|^ simplifies to 


c^7(n, '^2, a) K, c ri ■ r2 ■, 


(14) 


and we can explicitly calculate the kx distribution ([l^ , 
d^a^{q^qY) Oem C 


d(lnQ!) (PkT 

{q qY 

(i(lna) d'^k'T 


O-er). 

{kl + e^Y 

4 aem C a^(l — k'^ 


vr- 


{k^ + 


{217150*4+ [l + (l-a )*](4 + e*)} , ( 15 ) 

( 16 ) 


+2*14 


Note that for small a (p!5[) and (^) vanish like a^. This could have been expected since 
electomagnetic bremsstrahlung is known to be located predominantly in the fragmentation 
regions of colliding particles rather than at midrapidity. 

In the case of a nuclear target the transverse momentum distribution has to be modihed 
by eikonalization of ([I^) (see Appendix A), 

d?a^{q gq) 1 


(i(lna) d^kx (27r 


- / dVi (iVjexpliiiT)?! - fj)] ^Jo.f,) ^^(a.fj) E.,(fi,f2,Q) , 


( 17 ) 













where 


S^(fi,r2,a) = I d^b\l + exp^--aqq[a{n - f2)]T{b) 


— exp 


~aqq{ari)T{b) 


— exp 


~aqq{ar2)T{b) 


(18) 


The fluctuation wave functions in ([T7|) can be represented using (j^) in the form 


H a) a) = 

inj 




2 7r^ 


|m^ Ko{eri) Ko{er2) 


1 + (1 — a)^ 




in, f 


^ vr?* (1 - af Ko{eri) i^o(er2) , 

TT^ ' 


rir 2 
ri r 2 

\2 


A',(er,)A'i(er2)} , (19) 


( 20 ) 


where we average over the initial quark polarization and sum over the final polarizations 
of quark and photon. 

At first glance, one could think that the distribution is not modified by the nucleus 
in the case t/ 3> Ra, since the fluctuation is formed long before the nucleus and the quark 
interact. This is, however, not the case. Due to color Altering |jT^ the mean size of qq 
dipoles surviving propagation through the nucleus decreases with A. Correspondingly, 
the transverse momentum of the photon increases. In other words, a heavier nucleus 
provides a larger momentum transfer to the quark, hence it is able to break up smaller 
size fluctuations and release photons with larger kr- 

Note that one can also calculate the distribution with respect to the transverse mo¬ 
mentum pt of the final quark integrating the differential cross section over the photon 
momentum kr- The result turns out to be the same as ([T^) and (|T7|) with the replacement 
CK 1 — CK. 

We also calculated the nuclear dependence of the differential cross section (O- (0) 
using the dipole approximation for aqq{r). The details of the necessary integration can be 
found in Appendix B. As usual, we approximate the cross section by an A"'-dependence. 
The power n is then defined by 

d |ln (Ra^{q ^ q'y)/ {d{\ na) } 


n{kT, a) = 


d(lnA) 


( 21 ) 
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This power can also be A dependent. We performed calcnlations for A = 200. To simplify 
these calcnlations, we nsed the constant density distribntion, pa(^) = PoQ{Ra — r) with 
po = 0.16 fm“^. 

First of all, we calcnlated n{kT, a) for Drell-Yan lepton pair prodnction at 171 ^ = 
4 GeV. The results are shown in Fig. ^ for transversely and longitudinally polarized virtual 
photons (the two components can be extracted from the angular distribution of the lepton 
pairs). We see that n < 1 for < 1 GeV, i.e. the Drell-Yan pair production is shadowed 



Figure 3: The exponent (|2^ of the atomic number dependence parame¬ 
terized as A"' versus kx and a for transversely (left figure) and longitu¬ 
dinally (right figure) polarized virtual photons. 


by the nucleus. The shadowing is stronger for larger a 0 - Shadowing in the Drell-Yan 


process was first observed by the E772 Gollaboration |^. Their effect is, however, much 
weaker which can easily be explained because for Fermilab energies the radiation time (|I|) 
is quite short compared to the nuclear radius. This fact is taken into account in by 
means of nuclear formfactor. Then the data can be described quite nicely. (See also [^ .) 

An interesting result contained in Fig. is the appearance of an antishadowing region 
for kx > 1 GeV. This is the first case in which the coherence effects enhances rather than 
suppresses the radiation spectrum. It originates from an interference effect which is not 
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noticeable for the integrated qnantities. 

Nnclear antishadowing is especially strong for longitndinal photons and ~ 1-5 — 2 
GeV. Color hltering in nnclear matter changes the angular distribution of Drell-Yan pairs 
and enhances the yield of longitudinally polarized dileptons. The nontrivial behaviour 
of n for longitudinal photons at small kT is due to the dip at = 0 in the differential 
cross section for a nucleon, see Eq. (PIB]) . This minimum is hlled by multiple scattering of 
the quark in the nucleus leading to an increase of n{kT = 0) and a strong A-dependence 
of n{kT = 0). (Formally, for longitudinal photons n{kT = 0) goes to inhnity for A = 1, 
because the proton cross-section at kT = 0 is zero). 

Note that nuclear enhancement of Drell-Yan pair production at large kT was also 
observed experimentally [0. However, as was mentioned, these data were taken in the 
kinematical region of the Bethe-Heitler regime, i.e. t/ -C Ra- Therefore, they cannot be 
compared with our calculations. In fact the observation was explained quite satisfactory 
in [p. 


The /ct- dependence of n is expected to be nearly the same for different dilepton masses, 
down to the mass range probed in the CERES experiment at SPS CERN. However, the 
nuclear effects turn out to be quite different for real photons. Our results are shown in 
Fig. In order to compare with experimental dilepton cross sections and prompt photon 
production rates our results have to be convoluted with the quark distribution function for 
the projectile proton. Since the electromagnetic radiation steeply falls off with decreasing 
a (proportional to o?, see ([T^ ) - (|^)), the convolution effectively picks out large values 
of a where the nuclear effects are in turn expected to be large. Detailed calculations and 
comparisons with data are postponed to a later publication. 
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Figure 4: The same as in Fig. ^ but for real photons. 


3.2 Gluon radiation 


Now we can discuss bremsstrahlung in the non-Abelian case. Summing up the diagrams 
in Fig. 0 we get in impact parameter representation 
d?a^{q qg) 1 


d{hia) (PkT (27r) 


where (see Appendix A) 


2 . dVi dVs exp[i/cT(n - ^ 2 )] f) ^ 2 ) ^ 2 , of) , 


( 22 ) 


cr, 


,(ri,r2,a) = -{(Jgqq{fi,ri-ar2) + (Tgqq{f2,r2-ari)-aqq[a{fi-r2)]-agg{fi-r2)"^ . (23) 


Here Uggir) = | Uqqir) is the total cross section of a colorless gg dipole with a nucleon. 

Note that (^) reproduces several simple limiting cases: 

1. ) (Jg(ri,r 2 , a) vanishes if either of ri or r 2 goes to zero, which expresses the fact that 
a point-like quark-gluon fluctuation cannot be resolved by any interaction. To show this 
limiting behaviour one simply has to insert e.g. for F 2 = 0 the two relations o'gqq{ri, ri) = 
(Tggifi) and CTgqqiO, —ari) = <Jqq{—ari) = aqq{ari). (Quark and antiquark at the same 
point in space act like a gluon etc.) 

2. ) For a —> 1 the quark-gluon separation tends to zero and (p3D transforms into ([T3|) . On 
the other hand, at a —>■ 0 the quark-antiquark separation vanishes and (E^ takes again 


12 









the same form as (p!3D, except that the qq pair is replaced by a gluon-gluon dipole. 


If 19 

t^ 5 (ri,r 2 ,«)L«i = + (Tggir2) - agg[{n - f2)]\ = -a^.(fi,f2,a) 


Q:=l 


(24) 


We use the dipole approximation agg^rx) ~ Cr^, which is well justified in this case 


since the mean transverse quark-gluon separation is small at small a. In this case (^) 
and (ED lead to 


(Jg{ri,f2,a) 


+ 7(1 — ck ) 


C f 1 ■ f 2 


(25) 


This expression coincides with (0 up to the factor [1 -|- 9(1 — Q()/(4a^)]. Therefore, we 
can use the results (1^) - ([T6|) obtained for photon bremsstrahlnng which for ct ^ 0 lead 
to 

{q qg) K + 


(i(lnQ;) d'^kx 
d^a^{q qg) 


ckCI 


QC as k'^ + m^g 
TT^ {k'^ -I- 


(i(lna) d^kj 


12 C as rrdg k"^ 


(27) 


TT^ {k^ + rrdg)^ 

In contrast to photon bremsstrahlnng this cross sections do not vanish for a —»■ 0. This 
is a conseqnence of the non-Abelian nature of QCD []ro|. The radiating color current 
propagates throngh the whole rapidity interval between the projectile and the target 
providing a constant glnon density (^6]) - (|27|) with respect to rapidity. 


Eikonalization of the cross section ([2^) results in, 
d?a^{q —> qg) 1 


(i(lna) d'^kx 
where 


( 2 vr) 


dVi dVs e^v[ikx{ri - fa)] ^*gg{a, fi) ^gq{a, fa) Sg(fi, fa, a) , 


(28) 


Sg(fi,f2,a) = / exp[--ag-jQ((fi - fa)] 


-|- exp 


-^(^gg{ri-r2)T{b) 


— exp 


-2^gqq{ri,n-ar2)T{b) 


— exp 


-2^gqq{r2,r2-an)T{b) 


In the limit a -C 1, which is of practical interest at high energy 
form of (|2^) and Eq. (|29|) simplihes to 


(29) 

transforms to the 


Sg(fi,f 2 ,a)U«i = / d 6 <j 1 + exp[--agg(fi-f 2 )T( 6 ) 


— exp 


-\(^ggin)T{b) 


— exp 


~agg{f2)T{b) 


(30) 
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Note that the transverse momentum distribution for gluon radiation was calculated pre¬ 


viously in in the limit a —0 and rriq = nig = 0. Our results (|^) , (|30[) agree with 
that in in this limit. 


In (|30|) we make use of the fact that at zero qq separation a ^fgg-system interacts 
like a pair of gluons, o'gqq(f,f) = Cgg^r) = (9/4)(jgq(r). Therefore, (|^) - (|^) can be 
calculated in the same way as (0) - (0) in the electromagnetic case at a = 1 (see 
Appendix B), except that the fluctuation wave functions must be taken at a = 0. We 
assign an effective mass to the gluon, either of the order of the inverse conhnement radius, 
nig ^ 0.15 GeV, or in accordance with the results of lattice calculations for the range of 


gluon-gluon correlations of size nig = 0.75 GeV. We sum over the polarization of the 
emitted gluon. The numerical results are plotted in Fig. S. They are qualitatively similar 



Figure 5: The same as in Fig. but for gluons at a = 0 for different 
effective gluon masses. 


to those for photon radiation (see Fig. shadowing at small and antishadowing at large 
kx- However, the effect of antishadowing is more pronounced for light gluons. 

Antishadowing of gluons results in antishadowing for inclusive hadron production, 
which is well known as Gronin effect Although it was qualitatively understood that 
the source of this enhancement is multiple interaction of the partons in the nucleus, to 
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our knowledge no realistic calculation taking into account color screening was done so far. 
We expect that the Cronin effect disappears at very large kx, which would actually be in 
accordance with available data |2^. For a honest comparison with these data, one has to 
relate the kr of the gluon to that of the produced hadron, a step which lies not within 
the scope of this paper. 


4 Conclusions and discussion 


The main results of the paper are the following. 


The factorized light-cone approach [Q for the analysis of radiation cross sections 
was extended to treat the kx dependence, and was applied both to photon (real and 
virtual) and gluon bremsstrahlung. 


• The effects of coherence which are known to suppress radiation at long formation 
times, is only effective for small kx- At /ct > 1 GeV the interference instead actu¬ 
ally enhances the radiation spectrum. This was indeed observed for dilepton and 
inclusive hadron production off nuclei (Cronin effect). The enhancement of radia¬ 
tion by the coherence effects turns out to vanish at very large transverse momenta 
kr > 10 GeV. This was also observed in hadropro duct ion . 


• suppression and enhancement of radiation by the effects of coherence are quite 
different for transversely and longitudinally polarized photons. Both contributions 
can be separated by measuring the angular distribution of the produced dileptons. 


Note that we use Born graphs shown in Figs. S-l to derive expressions (|^) and others 
having a factorized form. As a result of Born approximation the dipole cross section 


u, 


qq{p) is energy independent. It is well known that the higher order corrections lead 


to a cross section rising with energy. HERA data suggest that this energy dependence is 
correlated with the dipole size tt- Therefore, the parameter ^(s) can be parameterized 
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as 


/ S \ 

^7(5) = Co - 
Vso/ 


(31) 


where sq = 100 GeV^, Co ~ 3. The power A^tt) grows with decreasing r^. This depen¬ 
dence is extracted from an analysis of HERA data in ^5 


Our results obtained for the radiation by a quark interacting with a nucleus are easily 
adapted to proton-nucleus collisions by convolution with the quark distribution in the 
proton. 

We plan also to extend our analysis to relativistic heavy ion collisions. The condition 
we use, 7/ 3> Ra is poorly satished at present hxed target accelerators, but are well 
justihed at RHIC or LHC. Indeed, if s^n is the total NN collision energy squared, for a 
gluon (photon) radiated at central rapidity, 

3 knp 


a = 


tf = 


\/snn 

^/SNN 


(32) 


(33) 


tun kx 

We conclude that at RHIC or LHC energies a <C 1 and that gluons with a few GeV 
transverse momentum are radiated far away from the nucleus, i.e. t/ 7§> Ra- Thus our 
calculations should be directly applicable. 
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Appendix A 


In this section we illustrate how to eikonalize the differential cross section in the case of a 
nuclear target and for the example of electomagnetic bremsstrahlung of an electron. The 
latter is described as propagating in a stationary field U{x), where x is a three-dimensional 
vector. 


The differential cross section reads, 

d^a 


ftp 


\M,i 


(A.1) 


dihia) d^pT d'^kx ( 2 vr)^ 
where kx and px are the transverse momenta of the photon and the electron in the final 
state. 

The radiation amplitude for a transversely polarized massive photon 7 * (o;^ = k‘^+m‘^,) 
has the form, 

Mjj = J d^x'^~\x,p 2 ) d ■ e* T^(x, pi) , (A. 2 ) 

where a = 707 are the Dirac matrices, and the wave functions ^(x, ^ 1 ^ 2 ) of the initial and 
final electron, are solutions of the Dirac equation in the external potential U{x), 

€ 1^2 — U{x) — m (3 + idV T(x,pi, 2 ) = 0 . (A.3) 

The upper indices ” —” and ” - 1 -” in (|A. 2 |) indicate that for the initial and final states the 
solutions contain in addition to the plane wave also an outgoing and incoming spherical 
wave respectively. 

If the energy is sufficiently high, 61^2 3> m, D it is natural to search for a solution of 
( |A.3| ) in the form of a polynomial expansion over powers of 1/e (e = 61 ^ 2 ); 

CX) 

, 

n —0 

T„(x,p) ~ e"'" . (A.4) 


Note that in the case of radiation of a longitudinally polarized photon it is sufficient 
to take into account only the main (Tq) which has a form, 

u{p) 


To(x,p)=e*^>^/(x,p/ 




(A.5) 
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where m(pi, 2 ) is the 4-component spinor corresponding to a free electron with momentum 
Pi, 2 , and the scalar function f{x,p) is a solution of the equation, 


(^A -I- 2ipV — 2eU{x)^ f{x,p) = 0 . 


(A.6) 


In the case of radiation of transversely polarized photons it is known that the two 
hrst terms in expansion ( |A.4|) are important. Their sum can be represented in the form, 


of Furry approximation, ^ 


To + = Ti. = I 1 - ^ \ /(f,^ 


(A.7) 


2e / ^ ■ 

One can estimate the accuracy of the Furry approximation using the following rela¬ 
tions, 

T - Tf = 5T = <F(T, ^ , (A.8) 

where <h(x, p) satishes the equation, 

IS. + 2ip ■ V — 2eU {x)-\-id ■ V U {x) =—^ {a-V U{x)^{d-V f{x,p)^ (A.9) 

It turns out that this correction to the Furry approximation for the electron wave 
function is of the order of f//e in the bremsstrahlung cross section. 

It is convenient (see below) to chose the axis along the momentum of the radiated 
photon. In this case one can represent the Furry approximation ([A.7|) for the functions 
T+(a;, pi) and T“(a;, P 2 ) in the form, 

m(Pi) 


T+(T,pi) = F+(x,pi) 

T]^(x,p2) = e*^^^4F-(f,p2) 


a /2 

U{P2) 

a/2 62 


where 


Di2 — 1 — i 


. a-V a(pi,2 - 


2 61,2 


2 61,2 


(A.IO) 

(A.ll) 

(A.12) 


k 


n = 


k ’ 

z = n ■ X 


Pi ,2 = \Pl,2\ ■ 
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In this case the functions F{x,p) and f{x,'p) are related as, 


F{x,p) = exp{ipx — ipz) f{x, p) 


(A.13) 


Therefore, F{x,p) = F^{x,p) has to satisfy the equation, 

(A + 2ip^ — 2eU{x)j F{x,p) . 


(A.14) 


The characteristic longitudinal distances in the problem under consideration xl ~ eirn? 
are much longer than the typical transverse distances xt ~ 1/m [^. Therefore, in the 
Laplacian A = (P/dz^ + {d/dx)"^ one can drop the hrst term d/^/dz^. Then (|A.14|) takes 
the form of the two-dimensional Schrodinger equation. 


d 




At 


F{x,p) , 

where p = \p\- We dehne F^ in accordance with the asymptotic behavior. 


(A.15) 


F+(x,pi) 

F~{x,P2) 


Z^Z—= — 00 


^ipirr 


^ip2Tr 


(A.16) 

(A.17) 


2:^2q_= + CXD 

Here we introduced new notations for transverse, r = xt, and longitudinal, z = xl, 
coordinates. 


It follows from (|A.15|) - (|A.17|) that these functions can be represented in the form. 


F~^{x,pi) = / d^riG{z,F, Z-,ri\pi) e^Frri ^ 


(A.18) 


F *{X,P2) = J d^r2G{z+,f2;z,f\p2)e ^ (A.19) 

where G(^ 2 , fb; A, a|p) is the retarded Green function corresponding to Eq. (A-15|) , 

G{z 2 ,f 2 -, zi,ri\p) =i5{z2 - zi) 5{f2 - n) (A.20) 


d A 2 / \ 


and satisfying the conditions. 


G{z2,r2] zi,fi\p) 

G{z2,r2;zi,fi\i^ 


Zi=Z2 
= 0 . 


= 6{f2-ri) 


Z1>Z2 


(A.21) 
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It is convenient to chose the axis 2 ; along the momentum of the radiated photon. Then 


Pit 

P2T 


kx 

a 

Pt - 


1 


a 


a 


(A.22) 


where kx and px are the transverse components of the photon and final electron momenta 
relative to the direction of the initial electron; a is the fraction of the light-cone momentum 
of the initial electron carried by the photon. 

We arrive at the following expression for the radiation amplitude ( A.2 ), 

1 


I d-^r^^rdzeM-^ At p2)G(z+,r,;z,r\p,) 
exp{i qrr,in z)T G{z,r; z_,ri\pi) exp{ipix n) , (A.23) 


2 p{l-a) 


where 


a m„ 


KJL t I U TTl 

^ 2(1 - a)Eg + ^ 

and Eg, rUg are the energy and the mass of the projectile quark. In the approximation 
considered in this paper when the fluctuation time substantially exceeds the interaction 
time, qmin -C and can be neglected. 

The vertex function in (|A.23|) reads, 


f 




— a u*{p 2 ) D^a ■ e* Di u{pi) 
ima {n X a) ■ e* + a {a X Vt) • e* 


i (2 - a) Vt • 


Xi 


(A.25) 


The operator Vt = d/dr acts to the right. xi ,2 are the two-component spinors of the 
initial and final electrons. 

In the case of a composite target the potential has to be summed over the constituents. 


U{f,z) = Y,Uo{r-ri,z - Zi) (A.26) 

i 

and the bremsstrahlung cross section should be averaged over the positions (?y, Zi) of the 
scattering centres. 
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Ml 


The averaged matrix element squared takes the form, 

OO OO 

= 2 Re y dzi J dz 2 J d^ri d^r[ d^r 2 d^r '2 d‘^r d?r' d'^p d?p' 

— OO Z\ 

X exp \ip2T {r 2 - ^ 2 ) - ipir W - n) - iqmin{z2 - Zi) 

X (^G{z+, ra; ^2, p\ P 2 ) G*{z+, r^; ^2, r' \ pa) 

X f'* (G{z 2 , p; zi,r\p 2 ) G*{z 2 , r'; zi, p' \ pi) 

X f {G{zi,f]Z-,fi\pi)G*{zi,p']Z-,T^^\pi)) , 


(A.27) 


where T' differs from T in (|A.25|) by the replacement 


^ d ^ , d 
V = — ^ V ' = ^. 
dr dr ' 

The following consideration is based on the representation of the Green function G in 
the form of a continuous integral p8[| , 

df{z) 


G{z2, ra; zi,fi\p) = J Vf{z) exp < ^ J dz 


dz 


i J dzU{f{z),z) \ , (A.28) 


where 


r{zi)=ri, r{z2)=f2, 


and the relation 
22 


J dzY, Uo{f{z) -fi,z- Zi) = Y x{r{zi) - fi) Q{z 2 - Zi) Q{zi - Zi) , (A.29) 


21 


where x(r^ = f)YdzUo(f, z). 

The mean value of the eikonal exponential is, 


^exp |i ^ [x {r{zi)) - x (f \zi))] ©(za - Zi)Q{zi - zi) | ^ = 


22 


exp<; --j dzn{z,b)a[f{z) -f\z)] ) , 


21 


where 


(A.30) 


cr(f-f') = 2 J d^p[l-exp{ix{r- ^ -ix{r' -^)] , (A.31) 
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and n{z^ h) is the density of scattering centres. 

Using these relations and performing integration by parts in ( |A.27|) , 
da^ 


OO OO 


dp 


d{\na)d‘^PTd'^kT (27r)‘^4p2 (^i _ a) 


- J dzi J dz2 j d^b d^pid^p2 


X exp 


— OO Z\ 

OO 


OO 2l 

iap 2 T p 2 - iapiT pi - J dzV{z, P 2 ) - J dzV{z,pi) 


Z2 


X 72 7i iu(^2,P2;u,pi Ip)- 


(A.32) 


The variables in this eqnation are related to those in ( |A.27| ) as, 

r 1 — r 1 


Pi = 
P2 = 


a 

ri - r2 

a 


b = -(Ui' + ri). 


Other variables in ( |A.27|) are integrated explicitly. 

Matrices 7 are related to T in (|A.25|) by replacement m am and d/dr^ d/dp. 


Absorptive potential V in (A-321) reads, 


V{z,^ = n{z,b)^{a ■ ^ , 
and W is the solntion of either of the eqnations, 

-^W{z 2 ,p 2 ]Zi,pi\p) = -^1 U(z2 ,P2;^1 ,Pi Ip) - U(p 2 , 2 : 2 )hU( 2 ; 2 ,P 2 ;^l,Pl Ip) , 

OZ2 2 a (1 — ajp 

(A.33) 

--^W{z2,p2]Zi,Pi\p) = ^-^^^—-^W{z2,p2]Zi,Pi\p)-V{pi,Zi)W{z2,P2-,Zi,Pi\p) , 

(A.34) 

with the bonndary condition 

W{z2,p2-,Zi,Pi\p) 

Using these eqnations and the relation. 


= ^(P2 -Pi) 


(A.35) 


Z2=Zl 


A(^ - £2 Ko{e\p\) = -27 i5{^ 


(A.36) 
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simple but cumbersome calculations lead to a new form for Eq. (|A.32|) , 


da^ 


a 


d{\na) d'^PT d'^kx (2vr)^ 


^ OO 

Re J dz J d^b d‘^pid‘^p2 d^p 


X exp 


OO Z 

iap2T p2 - id Pit Pi - J ^P2) - J dz'V{z!, pi) 


X '^^t{P2 - -V {z, pi) - V (z, P2)] ^t(Pi - ^ 

OO OO 

— 2Re J dzi J dz 2 J d% d^pid^p 2 d^p'-^^d^p '2 

OO z\ 

OO ' 

iap2TP2 - id Pit Pi - J dzV{z,p2) - J dzV{z,pi) 

Z2 —OO 

X I'kA -P 2 ') \V(Z2,P2) - V(^2,P2)] Ip) 

- l"(zi,Pi')l I'tIpi - Pi')|, (A.37) 


X exp 


where 


^t(p) = f K„{ep). 


(A.38) 


In the ultrarelativistic limit (p —>■ cxd) we have 


y {z 2 , p 2 ] zi, pijco) = 5{p2 - pi) exp 


22 

- J dzV{z,p2) 


2 l 


(A.39) 


The integrations over z, Zi, Z 2 in ( |A.37]) can be performed analytically, and we arrive at 
the expression 


da^ 


a 


dihia) d?pT d'^kx (27r)'^ 

X ' 0 T(n)' 0 r(h' 2 )S^(r,ri,r 2 ,a) , 


d^Ti d^r2 d^r exp i a f{pT + kx) + i{ri — r2) kx 


(A.40) 


where 


S^(r, ri,f2,a) = S^a(r + ri)^+S^a(f-f2)^-S(Q;r) - S^Q;(r+ fi - r2)^ , (A.41) 


S(p) = J d'^b <1- exp 


(^ip) 


T{b) 


(A.42) 
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and 





























The derivation of the correspondent expressions for glnon bremsstrahlung is done 


analogously. We skip the details and present only the results. 


Sg(f,ri,f2,a) = Ei(r,ri,f2,a)+S2(f,ri,f2,a)-S3(f,fi,f2,a)-S4(r,fi,r2,a) , (A.43) 


where 




S,(f,ri,f 2 ,«) = Jd%f^l-exp 

9 r / x 1 1 

cri(r,fi,f2,a) = - [cr(^r + (1 - a)f2j + (T(ri)J --(T(r + afi) ; 


9 

(T2(r,ri,r2,a) = - 


a 


(f - (1 - a)f2) + (T(r2) --a{r + ar2); 


(T 3 (f,fi,f 2 , a) = a(ar) ; 


(A.44) 

(A.45) 

(A.46) 

(A.47) 


/ \ 9 / \ 9 r 

a4(f,fi,f2,a) = o-(^f- Q!(fi - f2)j +-cr(^f + (1 - a)(ri - f2)j +-[cr(ri) + a(r2) 


4 V ^ ^ "V 8 L 

— a(r’+ (1 — afi + ar 2 ) — cr(^r — (1 — af 2 + ari) 


(A.48) 


This expression simplifies and gets the form of if one integrates in ( |A.4U|) over 
transverse momentum of the quark. Note that the last cross section < 74 (r, fq, fq, a) is the 
total cross section for a colorless system of two gluons, quark land antiquark interacting 
with a nucleon (compare with (p^ ). Here rq and rq are the transverse separations inside 
the qg and qg pairs and r is the transverse distance between the centers of gravity of these 
pairs. 


Appendix B 


In order to calculate Eqs. (|^) - (|I^) in the dipole approximation aqq = C r‘^, we need to 
evaluate integrals of two types: 

1 


h = 


(27r) 


rid^r 2 exp ikT (fq — fq) 


X Ko{eri)Ko{er 2 ) exp |- ^ (^frl + hr^ - 2grir2) | ; (B.l) 


and 


h = 


(2ir) 


d^rid^r 2 exp ik^ (fq — rq) 
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X 


{nr 2) 

rir2 


Ki{eri)Ki{er 2 ) exp |- ^ (/r^ + hrl - 251 ^^ 2 )| . (B.2) 


Here we use the notation, 


T{h) = j (/r? + hrl “ “^hnn) 


(B.3) 


We use the integral representation for the modihed Bessel functions, which reads 

, -2^2 
exD < —t — 


CX) 

/ X I f at I e^r 

^0 {er) = - J — exp {-t- 


At 


(B.4) 


^ 2^2 


1 / \ 1 f dt f 

-A'.(£r) = jy-exp|-t 

0 '' 


(B.5) 


After substitution of (E3) and ra into ( |B.1|) and ([B.2|) and making use of the 
following obvious relations, 

1 


h = 


d^rid^r 2 exp |ifcr(hi — n) 


4(27r)2 

- ^ + cr2 - 26rir2)| 

1 f fcf. (a + c — 26) 


(ac — 62) 


exp 


(ac — 62) 


(B. 6 ) 


h = 


d‘^rid‘^r2 {nr2) exp |ifcr(hi — r2) 


16(27r)2 
- ^ (a+ erg - 26 rir 2 )I 


1 


6 kl^{a + c — 26) 


{ac — 62)2 _ 52^3 


X exp 


/c|.(a + c — 26) 
ac — 62 


one arrives at. 


^ f dt du , . ^ 

Ji = - -exp(-M -t) Is, 

J t u 


h = 


o f dtdu , , ^ 

^ j 72 ^ exp(-n - t) h; 


where 


£2 £2 

a= — + /, c = — + h, b = g. 
t u 


(B.7) 


(B.8) 


(B.9) 
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Thus, for the general case in addition to the integration over the impact parameter 
one has to evaluate numericaly a two-dimensional integral over dt and du. 

The situation is simplified in the case of photon bremsstrahlung, when integration 
for the three exponentials in (|^) correspond to the following values of the parameters, 
respectively, 

f = g = 0, h= 2ca^T{h)] 

h = ^ = 0, /= 2ca^T{b)-, (B.IO) 

f = h = g = 2ca^T{b). 

In this case Eqs. (|B.7|) and (|B.7|) are reduced to one-dimensional integrals. 
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